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Probing Pseudogap by Josephson Tunneling 
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We propose here an experiment aimed to determine whether there are superconducting pairing 
fluctuations in the pseudogap regime of the high-T c materials. In the experimental setup, two 
samples above T c are brought into contact at a single point and the differential AC conductivity in 
the presence of a constant applied bias voltage between the samples, V, should be measured. We 
argue the the pairing fluctuations will produce randomly fluctuating Josephson current with zero 
mean, however the current-current correlator will have a characteristic frequency given by Josephson 
frequency uj — 2eV/h. We predict that the differential AC conductivity should have a peak at the 
Josephson frequency with the width determined by the phase fluctuations time. 
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One of the long-standing puzzles of the high- 
temperature superconductivity is the nature of the so- 
called pseudogap regime. The pseudogap regime occurs 
in the wide range of temperatures above superconduct- 
ing transition temperature in the underdoped cuprate 
superconductors. It is characterized by the suppressed 
quasiparticle density of states in the vicinity of the Fermi 
level. The similarity of the density of states in the 
pseudogap regime and in the superconducting state lead 
many [Q to believe that that the pseudogap itself is of 
a superconducting origin. In this view, the long range 
superconducting order in the pseudogap regime is de- 
stroyed by phase fluctuations. However, locally, both 
electronic pairing and fluctuating regions of supercon- 
ductivity should persist. Therefore, in this picture, the 
superconducting phase transition is believed to be a su- 
perconducting phase-ordering transition. Moreover, re- 
cent experiments by Orenstein and collaborators claim 
that local superfluid density is present even above T c in 
Bi2212 materials @. 

Whether or not the pseudogap indeed has a super- 
conducting origin remains to be verified experimentally. 
Such standard experimental test as the vanishing resis- 
tivity, or the Meissner effect are bound to fail. Both these 
test require spatial and temporal stability of the super- 
conducting phase on the time scale of the experiment. On 
the other hand, the pseudogap state can at best have su- 
perconducting order parameter that varies both in space 
and in time. A successful test may be possible if the 
fluctuations could somehow be stabilized by the proxim- 
ity to a "real" superconductor j|. Another approach is 
to probe superconductivity locally in space on the time 
scales comparable or shorter than the characteristic time 
of the phase fluctuations. 

In this letter we propose the first experiment of this 
type, which is based on the AC Josephson effect. The 
main point we make is that in the presence of phase fluc- 
tuations Josephson current j(t) across the tunneling con- 
tact is a random time-dependent quantity. It has a dis- 
persion that is the current-current correlator (j(t)j(t')), 



which is related to the corrections to the conductivity 
across the junction. This current-current correlator "re- 
members" about its Josephson origin and has a scale 
set by Josephson frequency and phase fluctuation time. 
Therefore, the conductivity of a junction will have a cor- 
rection due to current fluctuations, 



Act cx 1/V, 



(1) 



where V is the applied voltage across the junction. The 
crucial new aspect of the proposed approach is that we 
will focus on the characteristic time scale of frequency 
fluctuations, assuming that tunneling occurs in a small 
region where the spatial dependence can be ignored. We 
will focus on the time dynamics of the phase fluctuations 
in our analysis. 



2eV/h 




FIG. 1. Consider two parts of a superconductor brought 
into a tunneling contact at a temperature T > T c , so that 
the phase of the order parameter is no longer stationary. Our 
approach is to consider the Josephson current (Eq. pi) as a 
randomly fluctuating quantity, with the three major sources 
of randomness: the amplitude noise, the frequency noise, and 
the phase noise. We also assume that the contact is small 
enough for the superconducting parameters to be spatially 
uniform in the vicinity of the contact. 



When two pieces of a superconductor are joined by 
a weak link, a superconducting current begins to flow 
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through the link in the absence of applied voltage be- 
tween the superconductors. The current is related to the 
difference of the phases 0i and 02 of the superconductors, 

3 =josin(0i -<fa), (2) 

with the parameter jo related to the coupling strength 
between the two superconductors. For a superconductor 
in an equilibrium, evolution of the phase, 

d(f)/dt = -2fi/H, (3) 

is determined [^) by the superconductor chemical poten- 
tial, /z. Therefore, the phase as a function of time is 
0(t) = —2fJ,t/H + 0o, where 0o is the phase at time t = 0. 
For two coupled superconductors the phase difference 
evolves as 

Mt) - Mt) = 2(/i 2 - m)t/K+ 0i (o) - 2 (o). (4) 

The difference of the chemical potentials equals the ap- 
plied voltage, ^2—^1 = eV. Hence, if V = 0, both 
the phase difference and the current given by Eq. (j^) re- 
main constant. However, in the presence of a bias, V, the 
phase difference grows linearly with time and the current 
oscillates according to 

j(t)=j {uj)sm(tJjt + A(t>), (5) 

with the frequency wj — 2eV/h and the initial phase 
A0 = 0i (0) — 02(0). The effect of generating an alter- 
nating current by applying a constant bias to a supercon- 
ducting tunnel junction is called AC Josephson effect. An 
important feature of this effect is that the frequency of 
the generated current is only a function of the applied 
bias voltage and is independent of the microscopic and 
macroscopic parameters of the system. The scale of the 
frequency is about 0.5 THz per 1 microvolt. The AC 
Josephson effect is routinely observed in the supercon- 
ducting regime, ft can be observed either directly as the 
micro-wave emission from the oscillating current, or in- 
directly as "Shapiro steps" || in the DC I — V curves 
measured in the presence of the oscillating bias voltage 
component. 

For AC Josephson effect to be observable in the pseu- 
dogap regime, the measurement has to be local both in 
space and time. Suppose that above the transition tem- 
perature, superconductor can be modeled as a collection 
of superconducting islands of a characteristic size L, in- 
side which phase fluctuates at a rate A. Then, if the size 
of a contact between two superconductors is less than L 
then the superconducting state is essentially uniform in 
the vicinity of the contact. If the applied bias voltage 
is such that the Josephson frequency is larger than A 
then the Josephson oscillations are faster than the phase 
fluctuations dynamics, and hence can be approximately 
modeled as 



j(t) =f(uj,A,L) S m(L>jt + A<(>(t)), (6) 

with the amplitude j* being renormalized by spatial 
and temporal fluctuations of the superconducting phase. 
In general, the Josephson frequency can also fluctu- 
ate around its average value due to voltage fluctuations 
which are particularly important for small samples. How- 
ever, we assume that these effects can be absorbed into 
the overall fluctuations of the phase, A0(f). It is impor- 
tant to note that both A and L are functions of temper- 
ature, with L and 1/A diverging at the superconducting 
transition. The parameter L is related to the phase gra- 
dient correlation function W considered by Franz and 
Millis [Q. Relationship between A and L is a subject of 
an active interest. In the vortex diffusion picture, where 
phase fluctuations are produced by moving vortices, L 
is a distance vortex travels during the time 1/A, namely 
L = y/D/A. Here D is the vortex diffusion constant. 
This corresponding dynamical critical exponent is z = 2. 
Alternatively, if the phase fluctuations are governed by 
fast ballistic dynamics, the relation between the param- 
eters L and A should be LA = v* , where v* is propa- 
gation speed for the ballistic modes. This corresponds 
to z = 1. Using different geometries in the experimental 
setup that we propose below may help to determine the 
relevant model. 

A possible experimental setup that can be used to per- 
form the measurement of the AC Josephson effect in the 
pseudogap regime is shown in the figure 0. The crucial 
aspect is that the point of contact between the supercon- 
ductors be as small as possible. If the size of the contact 
becomes larger than L then in addition to the temporal 
phase fluctuations a the point of contact one needs to 
include spatial fluctuations, which can lead to a signifi- 
cant suppression of the effect. Another desirable feature 
is that the two superconductors be only a few afr-planes 
thick. This is because the size of the superconducting 
"islands" is likely to be more extended in the planes, 
compared to across the planes. Hence, we believe that 
the effect we propose is more likely to be observed in the 
geometry of figure [l|, although c-axis tunneling may also 
yield similar results. Finally, using very thin samples 
reduces the transition temperature 0, thereby making 
the pseudogap regime accessible at lower temperatures, 
where the thermal fluctuations are reduced. 

There are several ways the oscillating super-current in 
the pseudogap regime can be detected. Here we consider 
two methods: 1) differential AC conductivity measure- 
ments in the presence of constant bias voltage, 2) detec- 
tion of electro-magnetic radiation generated by the oscil- 
lating Josephson current. Although there is no coherent 
Josephson current in the pseudogap regime, the junction 
is expected to have strong response to the perturbations 
acting at the frequencies near ujj. Such super-current 
is also expected to generate a radiation peak at the fre- 
quency u) j, with a width of the peak governed by the 
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phase fluctuations. 

To make our qualitative arguments more formal we 
have to assign a particular form to the phase fluctua- 
tions. Here we make an assumption that the the phase 
difference between the superconductors, A(p(t) follows a 
diffusion process, as shown in Fig. |2|, with a variance 



<(A0 f -A^) 2 ) = 2A|t-£'| 



(7) 



and the initial phase A(f>o distributed uniformly in the 
interval [0, 2tt}. The factor of 2 appears because for a 
weak tunnel junction the phases of on the both sides of 
the junction fluctuate independently, each at a rate A. 




FIG. 2. We assume here that the phase difference, Acj>t, 
between the superconductors forming the Josephson junction 
follows a one-dimensional geometrical Brownian walk as a 
function of time, t. Such process is defined by a linearly in- 
creasing with time dispersion, Eq. {]?]). The initial phase, 
A(j>Q, is also a random quantity evenly distributed in the in- 
terval [0, 27t]. 



Viewing the Josephson current of Eq. (g) as a random 
quantity with a zero mean, we can characterize it by its 
dispersion and autocorrelation. The autocorrelation, ac- 
cording to the Kubo formula, determines the correction 
to the conductivity due to the fluctuating Josephson tun- 
neling, 



l r* 

Aa(tu) = — \ 
uv J_, 



j(t-t') 



(j(*)i(*')X 



(8) 



where brackets correspond to the </>-averaging, and aver- 
aging over time t is implied. Volume v is necessary for 
normalization. Substituting expression for the current 
from Eq. (||), we obtain 



•*2 



Act = 



2uv 



^(*- t ')(cos(wj(f - t') + Afr - A(j> t ')- 



cos(uj j(t + if) + A(f) t + A$ t >))dif. 



Since A<p t + A<t> v = (A(j> t - A0 O ) + (Afo ~ A</> ) + 2A0 O , 
after averaging over A<fo in the interval [0, 27r], the sec- 
ond cosine disappears. To average over (A0 t — A<fi t >), we 
invoke the relation (exp(iu)} = exp(— (m 2 }/2), valid for 
any normally distributed variable u with a mean zero. 
Then after the integration we obtain 



Act = 



3' 



1 



A + ilu — iujj 



+ 



1 



(10) 



As expected, the the real part of the conductivity, 



Re(Acr) 



J* 2 A 
4ojv 



(w - ojj) 2 + A 2 (uj + oj.j) 2 + A 2 



(11) 



which corresponds to the in-phase response, has two 
peaks near ±wj. The divergence as to — > has no 
physical meaning, since no superconductivity related re- 
sponse is expected on the time scales larger than the 
characteristic phase fluctuation time. This translates 
into the condition u> > luj for validity of Eq. (pd]). 
The imaginary part of the conductivity in the vicinity 
of loj is about two times smaller than the real part, 
and hence can be neglected in the total conductivity 
Abs ct = v/(Re ct) 2 + (Im ct) 2 . 

Therefore, we predict that if the pseudogap regime is 
superconducting in origin there should be a peak in the 
differential AC conductivity at the frequency luj oc V, 
with the peak value that scales as shown in Eq. (flu: 



Act oc j* 2 /Aujj oc l/V. 



(12) 



(9) 



This is the main result of this paper. While this correc- 
tion may be small relative to the normal (single-electron) 
current component, it can be extracted from the back- 
ground conductivity due to its extremely high sensitivity 
an applied external magnetic field. As is evident, the 
magnitude of the correction is inversely proportional to 
the phase-breaking rate, A, and as a consequence should 
be more easily observable at temperatures close to the 
superconducting transition. Consequently, a possible ex- 
perimental approach is to start arbitrarily close to T c and 
to measure the microwave radiation from the weakly de- 
phased Josephson current, and or to measure the differ- 
ential AC conductivity as proposed above. Then, gradu- 
ally incrementing the temperature one can probe how the 
spatial and temporal fluctuations of the order parameter 
phase grow with the temperature. 

In fact, a similar fluctuational AC Josephson effect can 
be searched for even in the conventional, "low-T c ," super- 
conductors, in the so-called paraconductivity regime [ p^[ . 
The paraconductivity regime is characterized by super- 
conducting order parameter fluctuations above T c , and 
experimentally is associated with the rapidly decreasing 
(but finite) resistivity in the vicinity of T c . Using the 
experimental setup we propose here, one could attempt 
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to study the dynamics of the dephasing timescales in the 
close proximity of T c in the paraconductivity regime. The 
difference between the conventional paraconductivity ef- 
fect and the pseudogap is that the pseudogap is believed 
to extend far beyond the paraconductivity range where 
the rapid changes in the material resistivity occur. 

Let us examine now more closely the assumptions that 
lead to the expression for the Josephson current, Eq. (||). 
Within the standard theory fTofl , the Josephson current 
is 



j(t) = 2elm[e-^*$ ret (el/)], 



(13) 



where the retarded correlation function <& Ic t(eV) can be 
obtained from the Matsubara correlation function 



2VV-k.-P I" dT'F^k : T,r')F(p,r' 7 r), 

kp J r-P 



via analytical continuation iuj — > eV + iO. Here Tk iP 
is a matrix element for tunneling from a state k on 
one side of the junction into a state p on the other 
side of the junction, and F(k, r, r') = (T T c k j( T )C_kx(T')) 
is an anomalous time-ordered Green functions. In Eq. 
(|l3|), we do not include the regular single electron con- 
tribution, proportional to G(k, r, r')G(p, r', r). The rea- 
son is that it does not carry the superconducting phase 
information and, therefore, does not produce the res- 
onant features away from zero frequency. In the ab- 
sence of phase fluctuations F is only a function of the 
time difference, F°(k,r - r') = (n F (-E k )e- E *\ r - T '\ - 
ni?(£ , k)e £k ' r ~ T ')/2i?k. Phase fluctuations can be incor- 
porated phenomenologically into the anomalous Green 
functions as phase factors 



F(k,r,r')^F(k,r,r')e 



(14) 



The form of the function </>(r, r') depends on the model 
of the phase fluctuations. Here we assume that 



(r,r') = 



6'(r - r') 



(15) 



where 4>{t) and 4>'{t — r') are uncorrelated Brownian mo- 
tions. Similar statistical properties of -F(k, r, r') can be 
obtained |ll| from a gauge transformation of the elec- 
tron operators, Ct = fAte 20 ^, under which F(t,t') — 
(T T ^ {T) ^ l {T'))e^+^'\ Since 9(r) + 9(r') = 
20(t) — (O(r) — 0(t')), and under realistic assumptions 
(0(r) — 6(r')) is only a function of (r— r' Lthis approach 
yields an expression equivalent to Eq. (|15|), except for 
the correlations induced between the functions 4>{ T ) an< i 
4>'(t — t'). In what follows we assume for simplicity that 
the correlations are absent. Then doing the average over 
the Brownian random process </>' and integrating over r' , 
for the Josephson current we obtain 



which is identical to the form of the Josephson current 
conjectured in Eq. (||). The function jo(z) is the ana- 
lytical continuation of the function jo(wj) which deter- 
mines the amplitude of the Josephson current in the ab- 
sence of the phase fluctuations. In the case of s-wave 
superconductivity with a constant gap A, this function 
is jo(eV) = (<T A/e)K(eV/2A), defined in terms of com- 
plete elliptic integral K(x). In the case of d-wave su- 
perconductor jo(eV) is also a nontrivial function of the 
relative orientation between lattices in two crystals. Its 
specific form is not important for our discussion. Finally, 
we should mention that in the current-current correlator, 
both </> and <t>' averages should be done on the product 
of currents, while we have done the averaging over 4>' 
independently in j(t) and j(t'). The qualitative results 
for conductivity, however, remain the same with the two 
peaks at the frequencies ±u>j. 

In conclusion, we propose to test the relevance of the 
phase fluctuations scenario in the pseudogap regime of 
the high-T c superconductors by investigating fluctuating 
Josephson current at T > T c . We focus on the temporal 
fluctuations of the phase assuming small-contact tunnel- 
ing to ignore spatial dependence of the phase. We ar- 
gue that although phase fluctuations will yield zero mean 
Josephson current, its autocorrelation function will pro- 
duce finite correction to the conductivity of normal cur- 
rent across the junction. AC conductivity will exhibit the 
peak at Josephson frequency u>j = 2eV/ h with the width 
determined by the characteristic phase fluctuation rate 
A. possible experimental test could be to measure the 
junction AC conductivity u{lj) in the presence of con- 
stant bias V and determine if it has a peak at u>j. We 
predict specific dependence Sa(cuj) oc V~ x of the peak. 
Specific experimental set up is shown on Fig. 1. 
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